Draft version January 30, 2013 

Preprint typeset using LAT^X style cmulateapj v. 5/2/11 



MODELING THE NUCLEAR INFRARED SPECTRAL ENERGY DISTRIBUTION 
OF TYPE II ACTIVE GALACTIC NUCLEI 

Paulina Lira, Liza Videla 

Dcpartamento dc Astronorma, Universidad de Chile. 

Yanling Wu 

Infrared Processing and Analysis Center, California Institute of Technology. 

Almudena Alonso-Herrero 

Institute dc Ffsica de Catabria, CSIC-UC, 39005 Santander, Spain 

David M. Alexander, and Martin Ward 

Department of Physics, Durham University. 
Draft version January 30, 2013 

ABSTRACT 

We present results from model fitting to the Spectral Energy Distribution (SED) of a homogeneous 
sample of Seyfert II galaxies drawn from the 12/im Galaxy Sample. Imaging and nuclear flux mea- 
surements are presented in an accompanying paper (Videla et al., 2013). Here we add IRS Spitzer 
observations to further constrain the SEDs after careful subtraction of a starburst component. We 
use the library of CLUMPY torus models from Nenkova et al. (2008ab) and also test the two-phase 
models recently produced by Stalevski et al. (2012). We find that photometric and spectroscopic ob- 
servations in the mid-IR (A > 5/im) are crucial to properly constrain the best-fit torus models. About 
half of our sources show clear near-IR excess of their SEDs above the best fit models. This problem 
can be less severe when using the Stalevski et al. (2012) models. It is not clear what is the nature 
of this emission since best fitted black body temperatures are very high (~ 1700 — 2500 K) and the 
Type II classification of our sources would correspond to a small probability to peer directly into the 
hottest regions of the torus. Crucially, the derived torus parameters when using CLUMPY models 
are quite robust,, independently of whether the sources require an additional black body component 
or not. Our findings suggest that tori are characterized by A/"o > 5, a > 40, r < 25, Zi > 40 degrees, 
Y < 50 and A° a ~ 100 — 300, where No is the number of clouds in the equatorial plane of the torus, 
a is the characteristic opening angle of the cloud distribution, r is the opacity of a single cloud, ZA, is 
the line-of-sight orientation of the torus, Y is the ratio of the inner to the outer radii, and A x ° & is the 
total opacity along the line-of-sight. From these we can determine that typical torus sizes and masses 
of 0.1 — 5.0 pc and 1O 4_6 M . We find tentative evidence that those nuclei with a detected Hidden 
Broad Line Regions are characterized by lower levels of extinction than those without one. Finally, 
we find no correlation between the torus properties and the presence of circumnuclear or more global 
star-formation. 

Subject headings: galaxies, infrared, seyfert, torus 



1. INTRODUCTION 

Evidence for the presence of circumnuclear obscura- 
tion in Active Galactic Nuclei (AGN) is undeniable. 
Obscuration has been determined from the excess 
a bsorption in X-r ay , UV and optical wav el engths 
((ILawrence & Elvij ll982t IMass-Hesse et all 119931: 
Turner et all 119971 : iMalkan et al.l 119981: iRisaliti et al l 
1999); see also further references in Videla et al., 2012; 
hereafter Paper I). However, whether this obscuration 
has similar properties in all sources or presents a wide 
variety of properties that vary from source to source is 
not yet determined. 

From a statistical view-point there is strong evidence 
that the properties of the obscuring material change as 
a function of luminosity, and possibly, redshift (Barger 
et al., 2005; Hopkins, Richards & Hernquist, 2007; Gilli, 
Comastri & Hasinger, 2007). Still, because of the higher 



spatial resolution, only local intermediate-luminosity (< 
10 46 ergs/s) sources can be studied in enough detail to 
disentangle the dust emission coming from the region 
close to the active nuclei from that coming from the re- 
maining host galaxy in the crucial near-IR domain. Also, 
the number counts in distant samples are dominated by 
Seyfert luminosity-class sources, and therefore their local 
properties might be a good representation of the highcr- 
redshift counterparts. 

All AGN, whether obscured or not, are thought to har- 
bor a central engine that emits because of the release of 
gravitational potential energy in an accretion disk sur- 
rounding a supermassive Black Hole. The Unified Model, 
proposed more than 25 years ago, is still strongly advo- 
cated to claim that, in fact, any AGN might look like a 
Type I source (i.e., unobscured) or Type II source (i.e., 
obscured), depending on the orientation of the obscur- 
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ing material, shaped as a torus, with respect to our line 
of sight (Antonucci 1993). If one assumes that all tori 
share similar physical traits, then the Unified Model im- 
poses some strong constraints on the observables, such 
as an expected correlation between the torus orientation 
angle and the optical depth towards the central source, 
among others. More recently, a new interpretation has 
been proposed, where Type I and Type II objects are 
preferentially drawn from the ends of the distribution of 
torus covering factor (Elitzur 2012). 

The early simple picture of a homogeneous torus (e.g., 
Pier & Krolik, 1992,1993; Granato & Danese, 1994; Ef- 
stathiou & Rowan- Robinson, 1995), which implied some 
clear differences in the properties between Type I and 
Type II sources, has been superseded by more com- 
plex models, where a clumpy media, with an overall 
geometry that still resembles that of a torus, give a 
much more realistic representation of the obscuring re- 
gion (Nenkova et al., 2002; Hdnig et al., 2006; Nenkova 
et al, 2008ab; Stalevski et al. 2012). Many efforts have 
already been carried out to fit the Spectral Energy Dis- 
tributions (SEDs) of Type I and Type II AGN using 
clumpy distributions of the dusty medium (e.g., Nikutta 
et al. 2009, Mor et al. 2009, Ramos- Almeida et al. 2009; 
Alonso-Herrero et al. 2011 (hereafter AH11); Mullaney 
et al. 2011). 

These new prescriptions imply that the correlation be- 
tween the different observables should be regarded as 
statistical ones, where individual sources might not be 
a good representation of the median of a sample. In the 
same way, the classification as Type I or Type II, also 
becomes a probabilistic problem. Therefore, samples as 
large and homogenous as possible have to be gathered in 
order to draw significant results about the torus proper- 
ties. 

With this idea in mind we have gathered near and mid- 
IR observations for a sample of Seyfert II galaxies to de- 
termine and fit their SEDs and in this way derive phys- 
ical and structural parameters of their tori. This paper 
is organized as follows: Section 2 presents the imaging 
and spectroscopic observations and their treatment pre- 
vious to the fitting process; Section 3 presents the clumpy 
models and the fitting procedure; Section 4 presents the 
obtained results and discusses the best-fit parameters; 
Section 5 looks into the possible correlations between 
the results and other characteristics of our sources, such 
as hydrogen column determined from X-ray data or the 
presence of a Hidden Broad Line Region (HBLR) inferred 
from spectropolarimetric observations. 

2. OBSERVATIONS 

The sample included in this work is formed by all 
Seyfert II galaxies in the southern he misphere found i n 
the Extended 12 /xm Galaxy Sample (|Rush et al.lll993l) . 
comprising 48 Seyfert II galaxietQ. 

The most important characteristic of the 12 /im Galaxy 
Sample for this work is that it includes a fairly unbi- 
ased sample of nearby (z<0.07) active galaxies that have 
been observed in many spectral regimes. It is selected 

1 This research has made use of the NASA/IPAC Extragalactic 
Database (NED) which is operated by the Jet Propulsion Labo- 
ratory, California Institute of Technology, under contract with the 
National Aeronautics and Space Administration. 



in the mid- infrared (MIR), minimizing possible biases: 
it includes elliptical, lenticulars and spirals galaxies (al- 
lowing to avoid systematic errors in the decomposition 
process of the surface brightness profiles, presented in 
Paper I); it includes a wide range in galaxy inclinations 
and most importantly, it includes a wide range of obscu- 
ration properties of the nuclear source (allowing a more 
general test of the nuclear emission models and probing 
a range of hydrogen columns of ~ 10 22 — 10 25 cm -2 ). 

The Type classification of the objects was ob- 
tained from existing catalogs of active galaxie s 
(jVeron-Gettv fc VeronMl99l"l : IHewitt fc Burbidgel fl99l . 
Some objects have been re-classified as modern observa- 
tions have become available. In order to corroborate the 
classification of the targets 38 objects were observed with 
the RC spectrograph on Blanco Telescope in CTIO, in 
two runs in August 2007 and February 2008 (see below 
and Paper I). 

2.1. Infrared Imaging and Optical Spectroscopy 

Details on the imaging observations for our sample are 
found in Paper I. We presented the near and mid-IR 
SEDs for 40 Type II Seyfert gal axies drawn from the Ex- 
tended 12 /im Galaxy Sample (Rus h et al.lll993l ). six of 
which had been determined previously by Alonso-Herrero 
et al. (2003). A detailed account on the data reduction 
and nuclear flux determination are given in Paper I. 

Additionally, optical long slit data were obtained in 
order to determine the nature of the nuclear ionizing 
source and the host stellar population and are also pre- 
sented in Paper I. This analysis has shown that 3 objects 
from the original sample have HII nuclei (MCG +0-29- 
23, NGC6810 and Mrk897) and therefore have been dis- 
carded from any subsequent analysis in this work leaving 
a total sample of 37 Type II Seyfert nuclei. 

2.2. Spectroscopic Observations 
2.2.1. Spitzer Data 

Most of the galaxies in our sample have been observed 
using the Infrared Spectrograph (IRS) on board Spitzer. 
The IRS provides moderate resolution spectroscopy from 
5.2 to 38.0 /im. It is composed of four separate modules, 
with two modules providing R^60-120 spectral resolu- 
tion over 5.2-38.0 /im (the Short-Low or SL, from 5.2 to 
14.5/im; and the Long-Low or LL, from 14.0 to 38.0 /im). 

The wide wavelength coverage of IRS spectra can allow 
a good representation of the mid-IR continuum of Seyfert 
galaxies. In contrast, ~ 10/im MIR ground-based spec- 
troscopy is normally dominated by Polycyclic Aromatic 
Hydrocarbure (PAH) features and the 9.7/im silicate ab- 
sorption present in that regime, resulting in very little 
free continuum to be observed. 

Results using Spitzer IRS data for nearby Seyfert 
galaxies have already been published by several authors 
(Buchanan et al., 2006; Deo et al., 2007; Melendcz ct 
al., 2008; Wu et al., 2009; Thomp son et al. 2009; Tom- 
masin et al., 2010). In more detail. IWu et all (2009) pub- 
lished the spectra for 103 AGN from the 12/xm Galaxy 
Sample, of which 44 are in this sample (the remaining 
seven lack Spitzer /IRS SL spectra). Most of the 44 ob- 
servations correspond to Program ID 3269 and were ob- 
tained between 2004 and 2005, a good match for the tim- 
ing of our photometry obtained between 2002 and 2004 
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Fig. 1. — Two extreme examples of subtraction of the emission 
from star forming regions using the templates published by Smith 
et al. (2007). From top to bottom in each panel the curves represent 
the original Spitzer-IRS data obtained using a small extraction 
aperture (§2.2.1) followed by the subtracted spectra using the four 
templates for star forming emission. The top panel presents the 
data for MCG-3-58-7, with very a small star forming component. 
The bottom panel presents the data for NGC 5953, which shows 
strong PAH features even after the subtraction attempts. This 
object has a very strong circumnuclear star forming region which 
clearly is contaminating the Spitzer observations. 

(see Paper I). The exceptions are IC 5063 and NGC 4941 
(PID 30572, observed in 2006 and 2007, respectively), 
and F 05189-2524 (as part of the IRS calibration cam- 
paign). 

For those objects observed in mapping mode (all 
sources in the 3269 program), sky subtraction was per- 
formed by differentiating the on- and off-source observa- 
tions of the same order in each module. In order to iso- 
late the AGN emission from the stellar components, and 
since rather large apertures were used to obtain the 1- 
dimensional spectra presented by IWu et all (|2009l ) , here 
we have reprocessed the data in order to have the small- 
est possible aperture, 2x2 pixels, or 3.6 x 3.6 arcsec 2 on 
the sky for the SL module. For data obtained with the 
IRS staring mode (IC5063 NGC 4941 and F 0518-2524), 
the reduction was done in the following manner. Indi- 
vidual pointings to each nod position of the slit were co- 
added using median averaging. Then on and off source 
images were subtracted to remove the contribution of the 
sky emission. Spectra from the final 2-D images were 
extracted with a point source extraction mode, which 
scaled the extraction aperture with wavelength to recover 
the same fraction of the diffraction limited instrumen- 
tal PSF (for details see Wu et al., 2009). Therefore, no 
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Fig. 2. — Examples of star- formation corrected Spitzer-IRS spec- 
tra compared with high-resolution ground-based, ~ 8 — 13 firm T- 
ReCS data for NGC 5135 and NGC 7130, as published by Diaz- 
Santos et al. (2010). The Spitzer spectra are presented before and 
after the subtraction of the starburst template and are shown with 
dotted lines. The subtracted IRS data have been scaled to match 
the T-ReCS data at ~ 10/jm. 

"small" aperture spectra are available for these 3 sources. 
All the spectra were flux calibrated using the IRS stan- 
dard star a Lac, for which an accurate template was 
available. 

2.2.2. Star Formation in the Spitzer-IRS Spectra 

As the apertures used to extract the Spitzer-IRS are 
larger than the spatial resolution of our photometry and 
it is likely to include extended emission coming from 
star formation located close to the nucleus, the spectra 
needs to be modeled to subtract a star-forming compo- 
nent before including the spectra in the nuclear SED of 
the galaxies in our sample. The presence of polycyclic 
aromatic hydrocarbons (PAH), which dominate the MIR 
spectra of star forming regions, will allow to isolate the 
A GN emission . 

ISmith et al.l (|2007t ) used low-resolution Spitzer-IRS 
spectra of a sample of nearby galaxies spanning a large 
range in star formation properties to construct 4 tem- 
plates in order to account for the differences of the PAH 
feature strengths and continuum properties. These 4 
templates were used to subtract the star formation com- 
ponent from the MIR spectra of our galaxies. Each tem- 
plate was scaled independently ensuring an optimal sub- 
traction of the PAH features but without oversubtracting 
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the continuum (i.e., wit hout yielding negat ive fluxes). As 
can be seen in Fig. 18 of lSmith et"alT(|2007l) the templates 
differ the most for A > 20/im. This region was excluded 
during the model fitting of our SEDs because the lack 
of photometric data at A > 10/zm makes it impossible 
to determine which template describes best the contin- 
uum emission for each particular galaxy for wavelengths 
longer than 20 /jm. 

In Fig. Q] extreme examples of the subtraction pro- 
cess are shown. The resulting spectra for MCG -3-58-7 
have no PAH features, and the subtraction only modi- 
fies the continuum for A > 20/im. The second example, 
NGC 5953, shows strong PAH residuals indicating that 
some emission from star forming regions is still present in 
the subtracted spectra. As discussed in Paper I, this ob- 
ject has a strong off-nuclear starburst region which is con- 
taminating the Spitzer data. Still, for the best subtrac- 
tion a large fraction of the PAH emission for A > 10/im 
is gone, while the 7.7^m feature leaves strong residu- 
als. Also, again it can be seen that the main difference 
between the four different subtracted spectra starts at 
A > 20/im. Only one object in our sample, MCG -2-8-39, 
did not require any correction for the presence of a star 
forming component since no evidence for the presence of 
PAH emission was found in its spectrum. 

We would like to asses the efficiency of our starburst 
component subtraction procedure. Using high-resolution 
ground-based T-ReCS data, Diaz-Santos et al. (2010) 
have studied the characteristics of the nuclear MIR emis- 
sion in nearby AGN. Two objects presented by Diaz- 
Santos et al. (2010), NGC 5135 and NGC 7130, are also 
common to our sample. 

Fig. [2] compares our subtracted Spitzer-IRS observa- 
tions with those obtained from the ground. The Spitzer 
spectra have been scaled so that the subtracted spec- 
trum has a common density flux at ~ 10/im with the 
ground-based data. 

Fig.[2]shows that a good match is achieved between the 
Spitzer subtracted spectra and the T-REcS data, partic- 
ularly for the regions dominated by the continuum. PAH 
residuals can be seen for both objects at the position of 
the 7.7 and 8.6/iin features, even though the ~ 12 — 13/im 
region presents a very clean subtraction. We can ten- 
tatively conclude that whenever a good subtraction of 
the PAH features is achieved the resulting spectrum can 
provide a good representation of the emission from the 
active nucleus. In fact, to avoid being too sensitive to a 
poor subtraction, whenever PAH residuals are seen, we 
use only a few spectral windows of the Spitzer observa- 
tions to constrain our model fitting (§4.1). It has to be 
stressed, however, that PAH and dust emission are not 
the only components present in the MIR spectra of star- 
bursts. Silicate absorption can be seen towards starburst 
nuclei with a mean flux ratio of the local continuum to 
the base of the line of 0.2 ±0.1 (Brandl et al, 2006). 

In particular, NGC 1144, NGC 5135, NGC 5953, 
NGC 7130, NGC 7496 and NGC 7582, show a strong star- 
burst component before the template subtraction, and 
significant structure afterwards. For NGC 1144 a clean 
PAH subtraction is achieved, but because of the intrin- 
sic silicate absorption towards starburst nuclei just men- 
tioned it is impossible to asses whether the residual dip 
seen at ~ 10/im is due to the presence of a dusty torus 
or not. NGC 7582 shows some PAH residuals at the po- 



sition of the 7.7 /im PAH feature after the template sub- 
traction, and therefore, again it is not simple to asses 
the presence of the residual 10/im feature. However, re- 
sults from small aperture, ground-based spectroscopy of 
NGC 7582 shows that in this case the presence of the sil- 
icate absorption is related to the nuclear source (AH11). 

2.2.3. Other Spectroscopic Data 

NGC 1068, was observed with the ISO satellite in 2001, 
and its spectrum was published by Le Floc'h et al., 
(2001), already divided into AGN and starburst com- 
ponents. Our photometry for this galaxy comes from an 
earlier date (1998 and 1999), but IR variability of this nu- 
cleus was of the order of 0.1 magnitudes in the -ftT-band 
during that period (Taranova & Shenavrin, 2006) 

NGC 7172 was observed using the TRecS spectrograph 
mounted on the Gemini South telescope (Roche et al., 
2007). Therefore this spectrum presents a narrower spec- 
tral coverage than the space-born observations. Because 
of the high spatial resolution, this spectrum isolates the 
AGN emission from the active nucleus. 

3. MODELING THE NUCLEAR IR SEDS 

3.1. Clumpy Torus Model 

As was suggested in many previous works, a clumpy 
structure of the distribution of gas and dust in the torus 
is more realistic than a continuous one. The fundamen- 
tal difference between clumpy and continuous density 
distributions is that radiation can propagate freely be- 
tween different regions of a medium populated by opti- 
cally thick clouds when it is clumpy, implying that cold 
dust may exist near the nucleus, and dust directly il- 
luminated by the central source may exist far from it. 
The difficulties in modeling such an environment plus the 
time consuming calculations (and technical limitations) 
prevented the developing of such approaching until re- 
c ently. 

iNenkova et al.l (|2002ft were the first in studying a 
cl umpy distribution for the gas and dust in the torus. 
In INenkova et al.l (|2008a| ) they presented the general for- 
malism for handling this clumpy media. The resulting 
S EDs (hereafter, the CL UMPY models) were presented 
in INenkova et al.l (l2008bD, with a n erratum on their cal- 
culations in INenkova "eFaLl (f20T0Tl . They assume that all 
the clouds are identical and characterize each one by its 
size (which should be much smaller than the separation 
between clouds), its opacity and its spatial distribution 
(the angular distribution of the clouds). Also the differ- 
ent dust temperatures in the illuminated surface and the 
dark side of the cloud are accounted for. The size distri- 
bution of the dust grains is that described by MRN, com- 
posed by the standard Galactic mixture (53% of silicate 
and 47% of graphite). Scattering is taken into account, 
dominating at short wavelengths (A < 1 /iin). They ar- 
gue that the illuminating spectrum makes no difference 
at wavelengths dominated by dust emission, A > 2 — 3/im, 
but at shorter wavelengths the AGN scattered radiation 
dominates. 

Through the clumpy treatment INenkova et al.l (|2002f ) 
can naturally explain the rather low and diverse dust 
temperatures found close to the nucleus o f NGC 1068 
paffe et all 120041: iSchartmann et al.l 120051 ). Further- 
more, it is found that the X-ray attenuating column den- 
sity is widely scattered around the column density that 
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TABLE 1 

Parameter values for torus models. 



Par Values Description 



CLUMPY Models (as described in Nenkova et al. (2008b)) 



~t 5, 10, 20, 30, 40, 60, V-band optical depth ol individual clouds 
80, 100, 150 

A/o 1-15 in steps of 1 Number of clouds at torus equator 

a 15-70 in steps of 5 Torus half-opening angle (0 at equator) 

q 0-3 in steps of 0.5 Radial cloud distribution exponent 

Y 5, 10-100 in steps of 10 Torus thickness (= R ut/Rin) 

IA 0-90 in steps of 10 Viewing angle (0 at polar axis) 



2pC Models (as described in Stalevski et al. (2012)) 



T9.7 0.1, 1, 5, 10, 15, 20 Total 9.7/im optical depth at equator 

Ri„ 0.5 Torus inner radius in pc 

Rout 15 Torus outer radius in pc 

p 0, 1 Radial clump distribution exponent 

q 0, 2, 4, 6 Polar clump distribution exponent 

© 50 Torus half-opening angle (0 at equator) 

Inclination, i 0, 40-90 in steps of 10 Viewing angle (0 at polar axis) 

Filling factor! 0.15, 0.25 Frequency of clumps 

Contrast 100 Fraction of dust in clumps 

Clump sizef 0.15, 1.2 in pc 

Resolution! lxlxl, 8x8x8 number of grid cells per clump 



Note. — j available parameter combinations arc: filling factor of 0.15, size clump 
of 0.15 pc and resolution of lxlxl; or filling factor of 0.25, size clump of 1.2 pc and 
resolution of 8x8x8. 



characterizes the IR emission, because the IR flux is col- 
lected from the entire observed area (averaging over all 
the clouds), while the X-ray opacity is calculated from 
one particular line of sight. 

The parameters for the CLUMPY models are the op- 
tical opacity of each cloud (ry), the number of clouds 
through the torus equator (A/o), the angular and radial 
distribution of the clouds (<r, measured from the torus 
equator, and r q , respectively), Y = R ou t/Ri„, and the 
angle between the axis of symmetry of the system and 
the line of sight (Zi). The values for each parameter are 
shown in Table [T] 

Essentially, the characteristics of the emerging emis- 
sion depend on whether the IR photons originate on the 
side of a cloud that is directly illuminated by the AGN, 
or the region that is heated by the radiation emitted by 
other clouds, the dark side. The optical depth will deter- 
mine how different the SEDs of these individual clouds 
are. The number of clouds of the whole torus popu- 
lation and their geometrical distribution will determine 
how many clouds are in the shadows of other clouds, 
while our orientation with respect to the torus will dic- 
tate how much of directly illuminated or dark parts of 
clouds we see. 

3.2. Two-Phase Torus Models 

Very recently Stalevski et al. (2012) have proposed two- 
phase models where dust might not only be found in 
clumps or clouds, but also in a diffuse medium filling the 
space between the clumps (hereafter, the 2pC models). 
The diffuse medium is controled by two parameters: the 
filling factor and contrast (see Table [IJ . The filling fac- 
tor determines the statistical frequency of clumps, with 
a value > 0.25 yielding enough clumps to form an inter- 
connected, sponge-like structure. The contrast parame- 
ter determines what fraction of the dust is found in the 



clouds, with a value of 1 corresponding to a smooth dis- 
tribution and a value > 1000 having effectively all the 
dust in the clumps. 

Stalevski et al. (2012) claim that the largest differences 
when this diffuse medium is included can be seen as a 
larger flux output in the near-IR. 

3.3. Fitting Procedure 
3.3.1. x 2 Test 

We compare our observed SEDs with models where 
the torus is described as a clumpy dist r ibution of gas 
and dust as proposed bv iNenkova et ahl (2008a. b), and 
which also includes a diffuse component as proposed by 
Stalevski et al. (2012). 

In order to include the Spitzer data in the modeling we 
had to take into account the scaling of the spectra, the 
spectral windows to be used, and the associated errors. 

For those objects where we had a N-band flux mea- 
surement, the spectra was scaled to match the photo- 
metric value. For those objects without a measurement, 
the average of the scaling factors found for the objects 
with N-band photometry was used. This was found to be 
~ 0.06 ± 0.03, with the error corresponding to the stan- 
dard deviation of the determined scaling factors. Photo- 
metric upper limits are also shown in the SEDs but they 
were not used during the fitting process. 

To avoid using Spitzer spectra with a poor subtrac- 
tion of the starburst component, we masked out strong 
PAH residuals during our fitting procedure, Also, for 
very noisy spectral data, the windows were averaged into 
"photometric" points, assuming the corresponding error 
of the N-band measurement. For well subtracted, high 
signal-to-noise spectra, the whole spectrum was used 
during the model fitting. We assumed a 5% error for 
each spectroscopic point (Wu et al. 2009) of the Spitzer 
observations. 
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TABLE 2 
Compiled Data for our Sample 
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Note. — f: For alternative galaxy names see Paper I. f: Spectral Class and Balmer Decrements (BD) from Paper I; unreliable BD values 
are labeled with §: Compilation of HBLR taken from Shu ct al. (2007); *: Data from Bassani ct al. (1999), Brightman & Nandra (2008), 
Greenhill, Tilak & Madcjski (2008), Noguchi et al. (2010), Panessa & Bassani (2002), Sazonov ct al. (2007), Shu et al. (2007), Tran (2003), 
Winter ct al. (2010); only direct X-ray flux components are quoted; I indicates Compton Thick sources; 0: Data from Thean et al. (2000) 
and Gallimote et al. (2006) at 8.4 and 5 GHz, respectively; radio flux for IC5063 at 5 GHz was obtained from Gregory et al. (1994); [\: 
Rl = fs I fsGHZ (see text for details); ^>: ncar-IR slopes from Paper I; values of a < 1 at a 1 sigma level are flagged with a (. 



We programmed a simple \ 2 routine that calculates 
the scaling factor needed shift the clumpy torus models 
to match the observed SEDs. The best fits correspond 
to the smaller \ 2 values. Notice that given our adopted 
scaling for the Spitzer data, the determined y 2 do not 
necessarily represent the actual goodness of the fit but 
allows to determine which is the best possible model. 

Some fits obtained using the CLUMPY library, how- 
ever, are at odds with the expected results. For example, 
the canonical Seyfert II galaxy NGC 1068, which has a 
well constrained SED, including a large number of IR 
photometric observations and a well determined MIR 
spectrum, is best fitted by models presenting a small 
value of the torus angle, this is, models with face-on 
orientations. This is not supported by a large body of 
evidence on the geometry of the central region in this ob- 
ject. In fact, based on this evidence, AH11 restricted the 
torus angle to values corresponding to high inclination 
only. 

It can be concluded, that the shape of the SED alone 
might not be sufficient to distinguish between models in 



the very large CLUMPY library. Whether or not a simi- 
lar situation would be found when using the 2pC models 
is not clear due to the much smaller set of available mod- 
els. 

3.3.2. L bo1 Estimates: Breaking the Degeneracy 

In order to break the possible degeneracy between dif- 
ferent CLUMPY models that present similarly shaped 
SEDs in the near and mid-IR regime, we have introduced 
another restriction. Different physical parameters or ge- 
ometry of the CLUMPY models, scaled to the same flux 
level, will necessarily imply different bolometric luminosi- 
ties of the central AGN. Therefore, we can use the 'ob- 
served' bolometric luminosities to discriminate between 
competing models. This is similar to the approach fol- 
lowed by Mor et al. (2009) when modeling Spitzer obser- 
vations of 26 luminous quasars. 

Determining the bolometric luminosity of Type II 
sources is not straightforward. The two most promis- 
ing luminosity proxies are the [OIII] line emission aris- 
ing from the Narrow Line Region (NLR), and the 2-10 
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Fig. 3.— Fits to the Spectral Energy Distribution of NGC 1068 
with and without (dashed and continuous lines, respectively) a L 
prior (see §3.2.2 for details). 

keV emission in the X-rays. However, both methods 
present important drawbacks. [OIII] fluxes can be dif- 
ficult to correct for slit losses and aperture effects for 
very extended NLRs, while determining an accurate ex- 
tinction correction can be problematic when the Balmer 
decrement is affected by absorption line features from 
the underlying stellar population. Probably of less im- 
portance would be taking into account the fraction of 
UV flux intercepted by the putative torus, and deter- 
mining the contamination by a starburst component. 
At the same time, the 2-10 keV emission is strictly a 
proxy only for Compton-thin sources (i.e., those with 
Nh < 10 24 cm -2 ), and determining accurate unab- 
sorbed 2-10 keV fluxes can be difficult for highly ab- 
sorbed sources (Nh ~ 10 23-34 cm -2 ) with average qual- 
ity data. See the work by LaMassa et al. (2010) for more 
discussion on the subject of luminosity proxies. 

We have compiled extinction and absorption corrected 
[OIII] and 2-10 keV fluxes for our sample (Table C2). 
[OIII] measurements are available for all of our sources, 
and it can be seen that they generally agree within a 
factor of 3 or better. X-ray fluxes are only available for 
about 40% of the sample, mainly due to the fraction of 
Compton-Thick sources. 

Using the expressions derived by Marconi et al. (2004) 
we have estimated the bolometric luminosities from the 
2-10 keV X-ray observations (see Tabled]). For the [OIII] 
data we have followed the work of Lamastra et al. (2009) 
who compared extinction corrected [OIII] luminosities 
and 2-10 keV fluxes for a sample of Type II Seyfert galax- 
ies to find £2-iokcV ~ IOLoiii- Adopting the Marconi 
et al. (2004) bolometric correction for the X-ray data, 
Lamastra et al. (2009) then tabulated 3 different bolo- 
metric corrections for log(Loiii) in the ranges 38-40, 40- 
42, and 42-44 ergs/s. Instead, we have combined the 
[OIII] vs 2-10 keV correlation given above with the bolo- 
metric correction found in Marconi et al. (2004) to find 
the cubic expression: 

logiom/i© = 9-45 - 0.76£ + 0.012£ 2 - 0.0015£ 3 (1) 

with C = \ogL ho1 1 'Lq — 12. The results are reported 
in Table [2] where we give the average value of L hal when 
more than one [OIII] flux is available. From objects with 



multiple [OIII] observations we find a rms of 0.3 ± 0.05 
in fractional flux. This rms value will be used as the 
standard deviation of the [OIII] flux distribution in what 
follows. 

Table [5] shows that when both estimates of the bolo- 
metric luminosities are available, in about 50% of the 
cases the values agree within a factor of 5 or better. 
However, in the other 50% of the cases there are im- 
portant disagreements. It is significant that in all but 
two cases (NGC 1144 and NGC 7172) L^f u > L^ 1(SkeV , 
maybe due to an insufficient absorption correction of the 
X-ray flux since several of these sources have large ab- 
sorbing columns. Because of this, and the complete avail- 
ability of [OIII] fluxes for our sample, we will use -^oiii 
as a prior for the bolometric luminosity of our sources. 

From the models, the bolometric luminosity of the 
AGN illuminating the torus is found as Lj^odel = @ x 
iird 2 , where d is the distance to the galaxy, and O is the 
scaling factor needed to shift the model to the observed 
SED data points. 

We use £^,^1 to restrict the best fit models by apply- 
ing a Bayesian approach. We compute prob(model|SED) 
oc prob(SED|model) xprob(model|Lbol), where we as- 
sume prob(SED|model) oc exp — (x 2 /2) (i.e., the likeli- 
hood that the SED has been obtained from the model) 
and prob(model|Lbol) cc e Xp (-(£&% - L h * del ) 2 /2a 2 om ) 

as prior, with croni = 0.3 x ^oni- 

Figure [3] shows the case study for NGC 1068. When 
restricting iJ^ ol dcl to the value implied by the observed 
[OIII] fluxes, the correct picture emerges. When applied 
to the entire sample, we find that about ^20-25% of 
sources change significantly their best-fit model results 
when applying the L ho1 restriction. 

3.3.3. Fitting Sources with near-IR Excess 

In many of our sources a near-IR excess is observed. 
These SEDs turn upwards for A < 3/im, this is, they 
present an index a < 1 in the near-IR at a 1 sigma level 
(for A/a oc A q ). These objects are flagged in Table[3J No- 
tice that MCG +0-29-23, NGC 6810 and MRK 897, which 
are classified as star-forming nuclei and therefore are not 
included in the present analysis, also show a near-IR ex- 
cess (Paper I). Also, NGC 34 and NGC 5953 have a < 1 
(although NGC 34 is not at a 1 sigma level), but as ar- 
gued below (see also §4.1), it is very likely that these 
sources are contaminated by a starburst component and 
therefore are not flagged in Table [5] as presenting a 'gen- 
uine' near-IR excess. 

In Paper I several possibilities to explain the presence 
of the observed near-IR excess were proposed: contribu- 
tion from a compact nuclear starburst or a nuclear stellar 
cluster, emission from the accretion disk, emission from 
a compact jet, or emission from a very hot dust compo- 
nent that survives within the sublimation radius of the 
torus. 

The first scenario, the contribution from a com- 
pact starburst, was not supported by the analysis of 
diagnostic-diagrams from emission line ratios and the 
analysis of the observed stellar continuum (see Paper I). 
Strong PAH residuals in the Spitzer data not necessarily 
imply a strong nuclear starburst component because of 
the large apertures. However, the nuclear optical spectra 
presented in Paper I suggests that NGC 34 might have 
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Fig. 4. — Example nuclear infrared Spectral Energy Distributions 
and optical spectra from Paper I for 4 of our sources with near-IR 
excess. The best model fit to the infrared SED is shown together 
with the accretion disk component after extinction has been applied 
to prevent it from appearing in the optical regime (dashed and 
dash-dotted lines, respectively). The amount of extinction applied 
to the accretion disk is shown at the top left corner of each panel. 
Notice that only TOL 1238-364 was observed under photometric 
conditions (Paper I), implying that the absolute flux calibration of 
the optical spectra for NGC1144, IRAS0 1475-0740 and NGC7130 
might be underestimated. 



a strong nuclear starburst, while NGC5953 has a very 
strong circumnuclear star-forming region. 

A second scenario is the presence of a luminous nu- 
clear stellar cluster. While nuclear clusters are found 
in most galaxy nuclei and are characterized by a light 
weighted stellar age of ~ 10 8-9 years (Walcher et al., 
2006), their luminosities (absolute magnitude zab = — 13 
for the most luminous cases; Cote et al. 2006) are not suf- 
ficient to make a significant contribution to the observed 
near-IR fluxes in our sample. 

The last three scenarios would require the leaking of 
the radiation from the central region through the clumpy 
obscuring torus. Since the nature of a clumpy medium 
would allow for certain lines of sight to peer directly at 
the central engine, even for highly obscured sources, in 
what follows we will study these scenarios in more detail. 
As 12/31 sources show evidence for a near-IR excess, 
we can postulate that the clumpy structure leaves, on 
average, ~ 40% of the lines of sight free of absorption for 
Type II sources. 

The first working hypothesis is that the presence of 
a near-IR excess could be due to central disk emission 
piercing through the torus for the particular line-of-sight 
we have of the system. To test this hypothesis SED 
fitting was done using the whole library of CLUMPY 
models, which besides the pure torus emission also have 
'SED+AGN' models, in which a torus is combined with 
the emission from the central nucleus. The advantage of 
using these models to simply adding an arbitrary power- 
law component is that the normalization of the AGN 
emission, which in turn illuminates the obscuring torus, 
is treated consistently with that of the emission from the 
torus itself. 

A consistency check to the best fit 'SED+AGN' models 
is to test whether the accretion disk emission is expected 
to be seen in the optical, since our sample is character- 
ized by a strong optical stellar continuum and the lack of 
a power-law component. We find that only three sources 
require a moderate amount of redenning (Ay < 1) to pre- 
vent disk emission from showing in the optical observed 
spectra (see Fig. HJ. This extinction could correspond to 
lines of sight at a grazing angle to some clouds in the 
torus, or to extinction introduced at larger scales. In 
principle, therefore, it is possible to have a disk com- 
ponent appearing at near-IR wavelengths without a dis- 
agreement with the optical observed continuum. 

A second consistency check is to examine the observed 
hydrogen column densities. A direct view of the accre- 
tion disk should correspond to objects presenting small 
X-ray inferred hydrogen columns, since these two obser- 
vations would probe the line-of-sight to the central re- 
gion of the active nucleus. Of the 12 sources identified 
as having a near-IR excess, one lacks X-ray observations 
(CGCG 381-051). Of the remaining, 8 have fairly large 
absorbing columns (log(iVij) > 22.5 cm -2 ; see TableEJ, 
while IRAS 01475-0740, NGC3660 and NGC4501 have 
log(N H ) = 21.6,20.5 and 22.2 cm" 2 , respectively. 

This result could be explained in the context of a 
clumpy torus if the probability of the power-law emit- 
ting region being obscured would be lower by ~ 70% 
(8/11) than that of the X-ray emitting region. In fact, 
in the canonical picture of the central engine of an AGN, 
the X-ray emitting region probably does not extend fur- 
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ther than a few gravitational radii from the central Black 
Hole, while the optically emission coming from the ac- 
cretion disk should reside hundreds of gravitational radii 
away, if a classical Shaku ra-Sunyaev a-disk is adopted 
(jShakura fc SunvaevllT976l ). 

However, this argument should also work in the oppo- 
site sense: those sources with small values of Nh should 
show a high probability to also exhibit an accretion disk 
component in the optical and near-IR. From Table [2] 
we find that, besides IRAS 01475-0740, NGC3660 and 
NGC4501, 4 other galaxies which do not show a near- 
IR excess have N H < 22 cm" 2 : ES0 33-G2, NGC5506, 
NGC5995 and NGC7590. In summary, out of 7 sources 
with low Nh columns 3 have near-IR excess and 4 do 
not. Hence, there is a ~ 50% chance to have a direct 
view towards the accretion disk when there seems to be 
unimpaired access of the innermost region towards the 
central black hole. This result comes from small number 
statistics but still argues against our null hypothesis. 

Moreover, our sample is composed entirely by objects 
with Type II classification from the absence of broad 
emission lines in their optical spectra, this is, their BLR 
is to a large degree completely obscured to us in direct 
light. But from the analysis above, in ~ 40% of the 
sources we might have a direct view of the accretion disk. 
As before, this would suggest that the probability of see- 
ing a larger structure is smaller than that of seeing more 
compact regions. This strongly argues against the in- 
terpretation of the near-IR excess as emission from the 
accretion disk. It is still feasible that in some of our 
sources the accretion disk component is indeed seen in 
direct emission, but it seems very unlikely that this can 
explain the large fraction of sources with near-IR excess 
in our sample. Moreover, as it will be further discussed 
in §4.3, other works have already claimed the presence of 
an extra near-IR component in the SEDs of Type I AGN, 
where the accretion disk and the innermost torus is read- 
ily visible. This is, the component is required besides the 
disk and torus emission, as modeled by CLUMPY. This 
extra component has so far been accounted for using a 
black-body spectrum. 

Consequently we conducted a fitting process using the 
torus CLUMPY 'SED' models plus a black-body com- 
ponent with a free scaling parameter (as opposite to the 
previous 'SED+AGN' modeling where both components 
were jointly scaled), representing emission from a hot 
component such as hot carbonaceous dust grains surviv- 
ing within the sublimation radius of the silicate dust. 
We tested temperatures in the 1000-2500 K range using 
steps of 100 K. 

The models proposed by Stalevski et al. (2012) might 
offer a different solution. Hot dust located in the diffuse 
component might be found further out in the torus and 
therefore increase the chance of making a contribution 
to the near-IR even for Type II sources. Because of the 
much smaller library of 2pC models currently available, 
and to test whether the treatment of the diffuse compo- 
nent detailed in Stalevski et al. (2012) can account for 
the near-IR excess emission, we conducted a simple \ 2 
minimization without adding extra components. 

3.3.4. Fitting Radio-loud Sources 

We have compiled nuclear radio measurements for 
most of our sources from Thean et al. (2000) and Gal- 



limore et al. (2006) at 8.4 and 5 GHz, respectively (see 
Table [2]) . We computed a radio- loudness parameter as 
Rl = F 5 chz/Fb as defined by Kellermann et al. (1989). 
Thean et al. (2000) fluxes were taken to 5 GHz assum- 
ing S v oc ^ _0 2 , implying a correction factor of 1.1 for 
these measurements. Since we do not have a direct mea- 
surement of the nuclear B magnitudes, we have used the 
Lorn va hies found in Table [2] to determine B fluxes as- 
suming the relation defined by Marconi et al. (2004). 
Radio-loud AGN are defined as those with Rl > 30. We 
can see that 5 of our sources fulfill this criteria, while 
most are in the Radio-quiet regime (Rl << 10). 

As detailed in Paper I, 3/12 of the sources with near- 
IR excess are classified as radio-loud: IRAS 01475-0740, 
NGC4501, and NGC7496. We can try to estimate 
whether Synchrotron emission is responsible for the near- 
IR excess using a correlation between optical and ra- 
dio flux determined for unbeamed low-power radio-loud 
AGN derived by Chiaberge et al. (1999) after interpolat- 
ing to the J-band wavelength. 

The predicted fluxes for our sources are extremely 
small, between 5 and 7 orders of magnitude below the 
observed fluxes, and clearly will not explain the observed 
near-IR excess, unless strong beaming takes place. We 
cannot rule-out this last possibility. Radio variability 
would uncover beamed sources. 

We conducted the fitting process for the 3 candidate 
beamed radio sources using the torus CLUMPY 'SED' 
models plus a power-law with a free scaling parameter. 
The power-law would correspond to a beamed jet which 
can have strong near-IR emission during a 'high-state' 
(Bonning et al., 2009). 

We used the calculated near-IR slopes as a first guess 
and explored 10 values around a in steps of 0.1. 

3.3.5. Fitting Sources with Strong Silicate Absorption 

A few of our sources present significant silicate ab- 
sorption at 9.7 /Ltm while showing strong near-IR emis- 
sion, this is, showing a 'broad' SED. The most ex- 
treme silicate absorption features are seen in NGC 7172 
and NGC 7582, while NGC 1194, IRAS 04385-0828, 
NGC 4388 and NGC 5506 correspond to more moderate 
cases. 

Using Sio to quantity the depth of the silicate fea- 
ture (5*10 = In F\/F Ci \, where F\ corresponds to the flux 
at the deepest point of the absorption feature and F c> \ 
corresponds to the flux in the interpolated continuum) 

we find that NGC 7172 and NGC 7582 have S w 2.0 

and ~ —1.4, respectively, while the other 4 sources listed 
above have 5io values between ~ —0.7 and —1. 

We find that it is not possible to properly fit 
NGC 5506, NGC 7172 and NGC 7582 with the set of cur- 
rent CLUMPY models because a combination of deep sil- 
icate absorption and SEDs with strong near-IR emission 
is not available. Even though NGC 5506 has a moderate 
absorption feature it has a particularly broad SED and 
therefore the obtained fits are poor. 

This is because CLUMPY models can only produce 
Sio > — 1 for intermediate optical depths, large number 
of clouds and edge-on lines of sight. At the same time, a 
large number of clouds yields models with a pronounced 
decrease in flux below 20/^m, as can be seen in Fig. 6 
of Nenkova et al. (2008b) , therefore producing 'narrow' 
SEDs. This issue is seen in the fits performed by AH11 
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Note. — Model results from the fitting procedure using CLUMPY models. For each parameter the average (mean), the mode (mod), and the median 
(med) with 1-sigma confidence limits (CL) of the probability distribution are given. Because of the discrete nature of the parameter space, for very narrow 
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Note. — Continuation from previous page, f Temperature of the black-body (BB) or index of the power-law (PL) secondary component, f The 
value of Ay along the line-of-sight are calculated following the equation found in text. Therefore no associated error estimates are give. 
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to NGC5506, NGC7172 and NGC7582 (notice that we 
share common photometry with AH11 for NGC 5506 and 
NGC7172). From their Fig. 6 it is clear how the fits to 
NGC 5506 and NGC 7582 become much steeper in the 
near-IR once the spectral information from the silicate 
absorption is taken into account. Only Circinus in Fig. 5 
of AH11 presents a combination of SED shape and sili- 
cate absorption strength compatible with the models. 

On the other hand, Stalevski et al. (2012) show that 
their 2pC models can have a value of S\o as low as 
—2.4. Crucially, their diffuse torus component con- 
tributes more significantly to the near-IR spectral range 
than the CLUMPY models, producing 'broader' SEDs. 
However, as noted by Stalevski et al. (2012) this effect 
is more substantial for face-on orientations, while we ex- 
pect that a large fraction of our sources prefer edge-on 
orientations. 

Levenson et al. (2007) have argued that a deep sili- 
cate absorption feature can be obtained if the emission 
source is deeply embedded in an optically and geomet- 
rically thick dusty structure. This is a reminiscence of 
the old continuous torus models, with the consequence 
that broad SEDs are not recovered. Alternatively, deep 
silicate absorption can be introduced by a thick screen of 
cold dust located further away from the central source, 
as proposed by AH11. Given the rather small expected 
torus sizes, it is possible that an external screen of dust 
will introduce the absorption feature along the line of 
sight towards us. However, as before, this screen will also 
reprocess the near-IR emission into longer wavelengths, 
and therefore the final SED will look narrower that the 
observations. 

We used DUSTY (Ivezic et al. 1999) to explore fur- 
ther this scenario using CLUMPY torus models as the 
incident spectrum on a slab of cold or hot dust. We 
found that a SED with strong near-IR emission is only 
recovered when the temperature of the slab is as high 
as ~ 1000K and for moderate optical depths. However, 
this combination of parameters essentially retains the sil- 
icate feature of the incident spectrum (ry ~ 20 implies 
T 9.7 ~ 1) while adding the emission from the screen of 
hot dust to the near-IR range. Essentially, this corre- 
sponds to a scenario where enough hot dust is added to 
the SED to increase its near-IR output without altering 
significantly the torus emission in the mid-IR. This can 
be regarded as an analogous model to the one proposed 
in §3.2.3 to explain the up-turn in the near-IR emis- 
sion of many of our sources. Notice that although this 
scenario might seem similar to that proposed Stalevski 
et al. (2012), our added black-body component is com- 
pletely arbitrary and has no physical connection to the 
torus structure. 

In summary, we finally conducted the fitting to sources 
with deep silicate absorption and a broad SED as de- 
scribed previously: allowing the presence of a black body 
component with temperatures in the 1000-2500 K range. 

4. RESULTS 
4.1. General Results 

Fig. [S] presents the fitting results to the SEDs of the 36 
sources found in Table 2. 

The 2pC model with the lowest x 2 is shown with a 
long-dashed line. No additional components were added 



to these models. 

The CLUMPY model with the highest probability is 
shown in Fig. [5] using a thick short-dashed line which 
includes the sum of any additional component, if present. 
The additional component (black-body or power-law) is 
separately shown using a thick dash-dotted line. 

In most cases, the probability distribution of the re- 
maining CLUMPY models tends to be highly peaked, 
in the sense that the probabilities rapidly fall to values 
much lower than that of the best-fit model. We have 
divided the probability distribution into three ranges: 
top-50%, 50 - 10% and 10 - 1% of the best fit prob- 
ability. The number of models found in each of these 
ranges is shown on the top-right corner of in each panel 
in Fig. [5j At most twenty models are shown in dark- 
gray, medium-gray and light-gray thin-continuous lines, 
for each of the probability ranges. If more than 20 mod- 
els are available for each range, then only 20 models are 
randomly drawn from the pool of models available and 
shown. This gives a representation of the level of depar- 
ture of the model SEDs from the best-fit. 

In what follows we discuss poor or inadequate fits. The 
mentioned sources will not be considered when drawing 
conclusions from the fitting process in §4.2. 

The need for mid-IR data to properly constrain the 
peak of the infrared emission is clearly seen in the fits to 
IRAS 00521-7054, ES0 541-IG12, and IRAS 01475-0740, 
where the best models only follow the rapid raise in near- 
IR flux but then flare-out for A > lOfim due to the lack of 
restriction in this range. These fits are also characterized 
by a much less peaked model probability distribution, as 
shown by the large number of models in the probability 
ranges defined above. Therefore all sources without mid- 
IR observations (this also includes IRAS 00198-7926 N 
and IRAS 00198-7926 S) will be excluded from the subse- 
quent analysis. Of these discarded sources, 2 have a < 1 
(see §3.3.3). 

NGC 34 and NGC 5953 show peculiar SEDs because of 
contamination from a starburst component, showing very 
strong near-IR emission. Both sources are also character- 
ized by Spitzer spectra showing the largest PAH residuals 
seen in the sample. We used a power-law to account for 
this component during the fitting process, but the results 
are poor. These two galaxies will not be considered in 
the analysis carried out below. 

NGC 4941 lacks enough data to constrain its SED and 
will not be analyzed any further. 

NGC 4501 presents a very steep SED throughout the 
near and mid-IR, possibly because of the presence of a 
beamed Synchrotron jet, as already discussed. This is in 
sharp contrast with the SED determined by Kishimoto et 
al. (2009a) using the Keck interferometer, which presents 
a canonical shape, rising from the near-IR to the mid-IR. 
In our observations, the presence of a dominant power- 
law component means that the torus properties cannot 
be successfully recovered from the fitting. 

In all, 27 sources will be considered in the following 
sections, unless otherwise noted. Of these, 14 required 
an extra component to account for the near-IR excess, 
i.e., very close to a 50% of the sample. More details are 
reported in §4.4. 

4.2. Results from SED Fitting & Best Fit Parameters 
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Fig. 6. — Top: Distribution of the best fitted values of q, Mo, o\ r, Zi, Y, logAv(los) and temperature of the additional black-body 
component if needed. The panels show the normalized histograms of the median (solid), and average and mode distributions (dashed lines) 
for the 26 objects with acceptable fits. Middle: Median distribution for objects that require a black-body component during their SED 
fitting (solid), compared with those that do not (dashed). Bottom: Median distribution for objects with a detected Hidden Broad Line 
Region (solid), compared with those without a HBLR (dashed). 



4.2.1. CLUMPY Fitting Results 

Table 3 presents the main results from the CLUMPY 
SED fitting procedure to the nuclear data of 27 sources, 
listing the weighted mean, the median, the mode and the 
67% confidence limits for the probability distribution of 
each torus parameter. 

Fig. [6] presents the distributions of the median (in 
solid), average and mode (both in dashed) parame- 
ters for the sample. Parameters are: q, A/o, c, r, Zi, Y, 
A v (along the line-of-sight) and black-body temperature 
when needed. In what follows we discuss the inferred 
values for each parameter. 

The exponent to the radial distribution of clouds shows 
a preference for q ~ 0. This corresponds to a flat distri- 
bution with the number of clouds showing a very weak 
dependence with the distance to the central Black Hole 
(Nenkova et al., 2008b). However, about half of the 
source require larger q values. 

The number of clouds along the equator shows a 
strongly rising distribution towards a large number of 
clouds, with A/o > 10 being clearly favored. This was 
also hinted by AH11. Nenkova et al. (2008b) showed 
that a rather small number of clouds is required by 
clumpy models to reproduce observed SEDs, with typ- 
ically Ao ~ 5. The difference with our findings could be 
due to our sample selection, which prefers obscured, IR- 
bright objects. Whether IR-bright AGN require a larger 
number of clouds in order to reproduce their properties 
is a tentative result from this work. 

A large number of clouds (Ao > 10) pushes the re- 
sults from the modeling of clumpy media to the limits of 



its parameter space (M. Elitzur, private communication). 
This is because the current calculations compute the ra- 
diation field produced by directly illuminated clouds and 
then solves for those clouds found in their shadow. The 
emission from clouds found in the shadow of other clouds 
is not taken into account. A complete solution would re- 
quire to iterate over the whole cloud population until 
a converging solution is attained. However, this is too 
demanding for current computer capabilities and these 
iterative steps are not followed (Nenkova et al. 2008a). 
Clearly, the larger the number of clouds, the larger the 
deviation between the current model prescriptions and a 
complete solution, particularly for compact geometries. 
For Ao — > oo the emission from a continuous torus should 
be recovered. 

The torus aperture angle shows a broad distribution 
with a peak at > 60 degrees and about 70% of the sources 
require a > 40 degrees. This is in agreement with AH 11 
who found that this is the case for Seyfert II galaxies, 
while Seyfert Is might show a narrower distributions of 
torus aperture angles. 

The optical depth of individual clouds in the V band 
is constrained to small values around r < 30, although 
some sources require larger values. 

The distribution of torus viewing angles shows that 
some objects are well fitted using intermediate values 
of Zi, but about half of the sample requires angles of 
70 to 90 degrees, reflecting the fact that our sample is 
IR-selected, and therefore should not be biased against 
heavily obscured objects. 

The distribution of torus thickness, Y, shows a ten- 
dency towards small values, with the bulk of the pop- 
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ulation requiring Y < 40. Nenkova et al. (2008b) ar- 
gued that given the level of isotropy observed in the IR 
emission of AGN, then torus should be rather compact 
(Y < 20) or present a steep (q ~ 2) rather than flat 
(q ~ 0) cloud distribution. We do not see this trend: 
sources with small and large values of Y all are more 
likely to require q ~ 0. 

The number of clouds along the line of sight can be 
determined using the expression: 

M((3) = Af -exp(~f3 2 /a 2 ) (2) 

where /3 is the angle between the torus equator and 
the line of sight (i.e., j3 = 90 — Zi). The product of the 
number of clouds and the optical depth of each cloud 
gives total optical depth of the torus along the line of 
sight in the V band: 

A v = 1.1- t-JV{/3) (3) 

Because of the rather narrow distribution of the values 
of r and jVq, the distribution of the total optical depth 
shows a strong peak at Ay ~ 30 — 300. 

The temperature of the black-body component is an 
additional parameter for the fits. The distribution of 
temperature values shows that most objects require a 
very high temperature (T > 1500K) and only 2 fits re- 
quire T ~ 1200 K. 

Crucially, the addition of a secondary component does 
not change the torus parameter distributions. Individu- 
ally, only a few sources show significant changes in some 
parameters after the black-body component is added (see 
Table 3). This can also be seen collectively when com- 
paring the histograms for objects where an additional 
component was required and for those where the torus 
model alone yielded a good fit (middle panels in Fig. [6]). 
The distributions are almost indistinguishable. 

4.3. Previous CLUMPY Fitting Results 

To check how robust our results are we compare 
our best-fit parameters with previous work which have 
studied common sourc es and which have also used the 
iNenkova etaLl (|2008afl CLUMPY models to represent 
the infrared SEDs of nearby AGN. Ramos-Almeida et 
al. (2009, hereafter RA09) and AH11 presented the mod- 
cling of Type I, Type II and intermediate Type Seyfert 
galaxies, of which 5 objects are in common with our sam- 
ple. AH11 adds N-band ground-based spectroscopy to 
the SEDs, while AR09 uses only photometric data. They 
both use a Bayesian inference algorithm (Bayesclumpy, 
Asensio-Ramos & Ramos-Almeida, 2009) to determine 
the best fit values to a given SED and the probabilis- 
tic distribution of the inferred model parameters. Notice 
also that AH 11 limited the range of Y to values up to 30 
for all objects, while the viewing angle was restricted to 
30 - 50 for NGC 5506 and 60 - 90 for NGC 1068, while 
our only restriction was the Lqiii P r i° r described in § 
3.2.2. 

For three of the objects presented by RA09 and AH11, 
NGC 1068, NGC 5506, and NGC 7172 the photometric 
measurements are common with this work (except for the 
Q-band measurements reported by RA09). IC5063 and 
NGC 7582, on the other hand, have independent pho- 
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Comparison of model parameters from R.A09, 

AH11, AND THIS WORK. 
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Note. — Notice that RA09 and AH11 limited the 
range of Y to values up to 30. AH 11 also incorporated 
a foreground host absorption component during their 
fitting corresponding to Ay = 7, 11, 8 — 13 magnitudes 
for IC 5063, NGC 5506 and NGC 7582, respectively. 



tometric observations. The ground-based spectra pre- 
sented by AH11 was obtained with much higher spa- 
tial resolution than the Spitzer and ISO data used here. 
However, examination of the spectra shows that the 
ground-base data are comparable to the space-born ob- 
servations. Hence the SEDs are not significantly modified 
by the different spectral data. 

Table @] presents the median and mode of the param- 
eter distributions determined by RA09 and AH11 and 
by this work for those objects common to the samples. 
Our approach gives very similar results to those found 
by AH11. The only significantly different model param- 
eters are q and A/q for NGC 1068. It can be seen that 
the largest differences are found between the work of 
RA09 and AH11, due to the inclusion of spectroscopic 
information around 10/im. This is clear indication of the 
importance of including detailed information of the SED 
around the silicate absorption feature. 

4.3.1. 2pC Fitting Results 
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results 


from the fitting 



procedure using 2pC two-phase models. As 
before, a ; indicates Compton Thick (CT) 
sources. 



Table [5] presents the results from the fitting using the 
2pC models for the 27 considered sources. The x 2 value 
from the best fit is also shown. About half of the sources 
have x 2 < 10 and we label these as acceptable fits. 

Examining Fig. [5] it becomes apparent that the 2pC 
models mostly fail to reproduce those SEDs that al- 
ready represented a challenge for the CLUMPY mod- 
els, namely, those with a near-IR excess and those with 
a strong silicate absorption accompanied by substantial 
near-IR emission. We have not attempted to obtain 
new fits using the current library of 2pC models adding 
a black-body component because of the rather limited 
range of parameter values available when compared with 
the CLUMPY models. 

Two parameters can be directly compared between the 
results obtained using 2pC and CLUMPY models: the 
total optical depth at the equator of the torus and the 
inclination angle subtended by the observer. Neither of 
these comparisons yield a proper correlation. In fact, 
2pC covers an optical depth range of T9.7 = 0.1 — 20, 
while CLUMPY considers values has high as T9.7 > 100. 

Fig. [7] presents the comparison between the best fit in- 
clinations angles. Even if only acceptable 2pC fits are 
taken into account, the distribution resembles a scatter 
plot. Again, notice that 2pC models do not cover inter- 
mediate angles (smaller than 40 degrees but larger than 
0) and the opening angle of the cloud distribution is fixed 
to 50 degrees. 

Therefore, it seems that the most important drawback 
from the 2pC models is the narrow range of parameters 
so far explored. 



100 



80 - 



a 60 - 
cv 

v 40 - 
tax) 



■ 

* ##* 

-x- -x- 



20 



20 40 60 80 
angle (CLUMPY) 



100 



Fig. 7. — Comparison between inclination angles obtained using 
the CLUMPY and 2pC models. All 27 sources are plotted with 
an asterisk. Sources with 2pC fits with values of x 2 < 10 are also 
shown with a circle. 

4.4. Sources with near-IR excess & strong Si absorption 

As already discussed, several sources were identified in 
Paper I as having a near-IR up-turn (i.e., a < 1; see also 
Table [2]) and their SEDs were fitted using an additional 
black-body component, which considerably improved the 
quality of the fits. The required temperatures are gener- 
ally very high (T ~ 1700 — 2500 K) and the emission typ- 
ically peaks at around 2/im. Only 2 sources require lower 
temperatures (T ~ 1200 K): NGC 4968 and NGC 5506. 

Several other works have found necessary to use an 
additional component to explain the excess of near-IR 
emission observed in luminous Type I QSOs with respect 
to the CLUMPY torus models (Mor, Netzer & Elitzur 
2009; Mor & Trakhtenbrot, 2011; Deo et al. 2011; Landt 
et al. 2011). All these works find that the black-body 
emission is characterized by T ~ 1200 — 1400K, while 
we find that the required temperatures are much higher, 
with T - 1700 - 2500 K. These values are too high to 
interprete this component as emission from hot dust and 
therefore its nature remains unexplained. Notice, how- 
ever, that given the difficulties in isolating the AGN near- 
IR emission in obscured Seyferts compared with the dom- 
inant AGN emission in luminous QSOs might contribute 
to the problem. 

Of the 12 sources with a < 1 in Table [2l 3 corre- 
spond to radio-loud sources and were also fit with an 
additional power-law component (§3.3.4). IRAS 01475- 
0740, however, had not enough mid-IR information to 
properly restrict the fitting. NGC 4501 seems dominated 
by a power-law and the inclusion of this additional com- 
ponent allows for a proper fit to the SED, as seen in 
Fig. [5] This component, however, heavily dilutes the 
torus emission, and it is not clear whether the derived 
parameters are representative of its intrinsic emission. 
The last radio loud source, NGC 7496 can be success- 
fully fit using a black-body or a power-law as secondary 
component. Fig. [5] shows the results from the fitting us- 
ing a black-body component, while Table 3 shows the 
best fit parameters for both cases. It can be seen that 
the results are consistent with each other. 

In §3.3.5 it was also discussed that sources presenting 
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deep silicate absorption were fitted using an additional 
black-body component to supply the near-IR flux lack- 
ing in CLUMPY torus models that have a strong 9.7/im 
absorption features, which is the case of NGC 5506, 
NGC7582 and NGC 7172. 

A good fit was obtained for NGC 5506 and NGC 7582, 
as can be seen in Fig. [SJ In Table 3 we report results for 
fits with and without the black-body component. For 
NGC 7172, on the other hand, no possible combination 
would reproduce both, the extremely deep 9.7/im ab- 
sorption feature and the near-IR photometry. However, 
the best fit 2pC model provides a reasonable fit for this 
source. It corresponds to a model with a 90 degree in- 
clination angle and a large T9.7 value. Large T9.7 corre- 
sponds in fact the regime where the two-phase medium 
introduces the largest departure from the clumps-only 
models. The CLUMPY results are also consistent with 
a large optical depth and high inclination angles. 

Fig. [5] shows that other sources require additional near- 
IR flux, other than those with a index a < 1 or strong Si 
absorption: NGC 4968, NGC 5995, and MCG -3-58-7. A 
very good fit was found for NGC 4968 with a CLUMPY 
model and a black-body component with a temperature 
of - 1200 K. The fits to MCG -3-34-64 and MCG -3-58-7 
are rather poor in the near-IR region, but the fit provided 
by the 2pC models to MCG -3-58-7 is quite good (x 2 = 
6). Table 3 presents the results using CLUMPY models 
with and without the additional black-body component 
for all these sources. 

It is very encouraging that the inclusion of a black- 
body component does not have a significant impact on 
the torus parameters derived from the fitting process us- 
ing CLUMPY models. This can be seen in the middle 
panel of Fig. [6j but also when examining individual fits 
reported in Table |U This is in contrast with the work 
of Deo et al. (2011), who found significant changes in 
the best fit parameters when introducing a black-body 
component to their fits. This could be due to the Type 
I nature of their sources and therefore a different level 
of constraining coming from the 9.7/im silicate feature in 
their sample of QSOs. The largest changes are seen in 
parameters q and Y. One of the most interesting findings 
from this work is that most QSOs require large inclina- 
tion angles and a small number of clouds. 

In summary, to solve for the lack of sufficient near-IR 
emission in CLUMPY torus models, we added an addi- 
tional black-body component to ~ 50% (14/27) of our 
sources and obtain acceptable fits in most cases. How- 
ever, the temperatures of these components are too high 
to correspond to dust emission, and therefore its true na- 
ture is unclear. 2pC models can provide better results for 
some of these sources, but the limited parameter space 
currently explored by these models does not allow to as- 
certain that the two-phase approach is a definite solution 
to the lack of near-IR emission in observed SEDs. 

5. ANALYSIS 

In the previous section it has been established that the 
results obtained using CLUMPY models are quite robust 
to the addition of a black-body component in the near- 
IR. Because of the large parameter space explored by 
these models and the well restricted results obtained for 
most parameters in §4.2.1 in what follows we will only 
consider CLUMPY results for our analysis. 



We have compiled ancillary data for our sample from 
the literature. Table [2] presents results on the detection 
of Hidden Broad Line Regions, [OIII] fluxes and lumi- 
nosities, hydrogen column densities determined from X- 
ray observations and hard (2-10 keV) fluxes, radio fluxes, 
and Balmer decrements as determined in Paper I. 

5.1. The Inferred Hydrogen Columns 

One of the most powerful diagnostics to characterize 
the different classes of AGN is the inferred hydrogen col- 
umn density, Njj, as determined from the photoelectric 
cut-off experimented by the power-law X-ray spectrum 
emitted by the central source. The observed values of Njj 
probes the amount of material along the line of sight to- 
wards the active nuclei and correlate strongly with other 
diagnostics to determine the Seyfert type. In fact, it has 
been shown that while Type I objects suffer from little 
absorption, Type II systems usually present absorbing 
columns of 10 22 cm -2 or more (Smith & Done 1996; Tur- 
net et al. 1997; Maiolino et al. 1998; Bassani et al. 1999). 

However, measurements of the column Njj towards the 
central region obtained through other methods, such as 
the ratio H a / Hp for broad emission lines (since narrow 
lines would probe the extinction affecting the much more 
extended narrow emission line region), typically give 
smaller values of Njj than those obtained from X-ray ob- 
servations (Maccacaro et al., 1982; Reichert et al., 1985). 
The optical depth can also be inferred using key features 
in the extinction curve, such as the 2200A 'bump' and 
the silicate absorption features in the mid-IR. 

Possible solutions for these differences have been pos- 
tulated: a dust free inner region (interior to the torus) 
could be responsible for the excess column probed by 
X-rays; an anomalous H a /Hp ratio could be due to the 
collisional effects present in the high-density clouds found 
in the broad line region; the line of sight probed by the 
X-rays could be significantly different to that probed by 
other estimators; Maiolino et al. (2001ab) and Gaskell 
et al. (2004) argue that AGN environments might have a 
different dust size distribution, either because of the pres- 
ence of larger grains or because small grains are depleted, 
although other works argue for normal dust properties in 
AGN (Mason et al, 2004; Nenekova et al., 2008b). X-ray 
absorption by dust-free material within the sublimation 
radius of the torus should not be significant as this ion- 
ized gas would produced very intense narrow emission 
lines which are not observed (Maiolino et al. 2001a). 

Shi et al. (2006) have shown that there is a broad but 
clear correlation between the strength of the 9.7/im sil- 
icate feature and the Njj columns derived from X-ray 
data. In very broad terms, unabsorbed systems show 
9.7/im feature in emission while absorbed systems show 
it in absorption. In Fig. [8] we compare the column den- 
sity Njj along the line of sight derived from our SED 
fitting with those obtained from the photoelectric cut-off 
from X-ray observations. Arrows show objects for which 
Njj > 10 24 ~ 25 cm~ 2 upper limits have been determined. 
One system has a Njj > 10 22 cm~ 2 upper limit coming 
from observations in the soft X-rays. No clear trend is 
found in the plot, and we are not able to reproduce the 
findings of Shi et al. (2006). Crucially, we do not find 
a systematic offset between the two measurements. We 
also show the combination of parameters No and r to 
yield the total Hydrogen column Njj (assuming a gas- 
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to-dust ratio N H /A V = 1.79 x 10 21 cur 2 mag -1 ) at the 
torus equator and along the line-of-sight. 

The presence of clumpy media around the active nu- 
cleus seems to be the best way to interpret our results. 
While our SEDs probe the average conditions of the 
dusty medium in emission, the X-rays probe a particu- 
lar line of sight towards the nucleus in absorption. Rapid 
and dramatic changes of the X-ray Njj in NGC 1365 and 
NGC4388 (Risaliti et al, 2005; Elvis et al., 2004) seem 
to validate this scenario (although Elvis et al. oppose it). 

So, can we consider our column densities more repre- 
sentative of the average, long term conditions of the phys- 
ical conditions of the duty torus? Probably yes, but these 
will have to be revised as new models and better obser- 
vations come along in the future. Since from the model 
fitting no Nh values are found outside the io 22 - 5-25 
cm~ 2 , this might be a reasonable range of columns to 
be adopted as representative of the average values for 
Type II Seyferts. 

Compton Thick nuclei are defined as those where the 
X-ray derived Nh columns are in excess of 10 24-25 cm~ 2 . 
The current compilation of Nh columns presented in Ta- 
ble [2] shows that of the total number of objects with 
measured column densities (30 out of the 39 galaxies in 
the full sample), 11 are Compton-Thick. This is in line 
with previous findings (Risaliti et al., 1999; Bassani et 
al., 1999; Bassani et al., 2006). 

Out of the 11 CT sources found in our sample, 8 have 
acceptable SED fits: the bona fide CT galaxy NGC 1068, 
NGC 1320, Tol 1238-364, NGC 4968, NGC 5135, F 15480- 
0344, NGC 7130 and NGC 7674. The inferred angle that 
the torus axis subtends with our line of sight are in the 
70-90 degree range for 7 objects, while for NGC 1320 the 
angle found is rather small (~ 20). See Table 3. 

AH11 has shown that the probability for photons to 
scape the obscuring material can still be low for sys- 
tems with intermediate inclinations, provided that the 
angle subtended by the torus is fairly large. This is the 
case for most of our Seyfert II systems which present 
rather large values of a. The exceptions are Tol 1238- 
364, IRAS 15480-0344, NGC 6890, NGC 7130, MCG-3- 
58-7 and NGC 7674. Interestingly, 4 of these sources are 
CT systems. It should be remembered, however, that 
the X-ray derived Nh columns measure a very specific 
line-of-sight property of the central region, as clearly val- 
idated by those objects with variable columns (Risaliti et 
al., 2005; Elvis et al., 2004). 

Elitzur (2012) has recently proposed that Type I and 
Type II nuclei are examples of objects preferentially 
drawn from the two ends of the distribution of torus 
covering factor, with Type II sources being examples of 
particularly puffed-up tori, characterized by large values 
of A/"o and a. This is clearly supported by our results. 

5.2. The presence of Hidden Broad Line Regions. 

The presence of a Hidden Broad Line Region is the 
most clear indication that at least some Type II AGN 
have broad emission lines. About 40% of Seyfert 2 galax- 
ies in the 12/im Sample are found to have HBLRs, in 
agreement with the fraction found in optically selected 
Seyfert 2 samples (Tran 2003). The fraction is closer to 
50% for the 27 objects with good SED fittings. Unfortu- 
nately, there are significant differences in the sensitivity 
of the spectropolarimetric studies found in the literature 
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Fig. 8. — Top panel: comparison of column densities along the 
line-of-sight, Njj (cm -2 ), derived from the CLUMPY emission 
models for the torus and derived through X-ray observations. Mid- 
dle panel: distribution of the number of clouds at the torus equator 
{Mo) versus the total Nfi column derived from our fittings, also 
at the equator. Continuous lines correspond to the total opacity 
for different optical depths of individual clouds (r = 10, 150, 500). 
Bottom panel: same as above but this time showing the number 
of clouds (Af(/3)) and the Nh column along the line-of-sigh. 
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and since the data are not provided in many of these 
works it is not possible to visually verify the presence or 
absence of the HBLRs. 

The reasons behind the lack of detection of a HBLR in 
some Type II sources, have been a matter of heated de- 
bate (Heisler, Lumsden & Baily, 1997; Alexander 2001; 
Gu et al., 2001; Tran 2001; Lumsden & Alexander, 2001; 
Gu & Wang; Tran 2003; Deluit 2004; Lumsden, Alexan- 
der & Hough, 2004; Zhang & Wang, 2006; Shu et al., 
2006; Haas et al., 2007). Some works advocate that the 
non-detection of a HBLR is due to the presence of a dom- 
inant galaxy component that dilutes the emission from 
the active nucleus. Others interpret the observations as 
evidence for the existence of a different class of AGN 
where a BLR is not present. 

It could be imagined that in objects that truly lack a 
broad line region, no dusty structure needs to be present 
either, but this does not have to be the case. The non- 
HBLR objects might have an active nucleus and a torus, 
lacking only the BLR. Evidence seems to support this: 
Haas et al., (2007) looked at the mid-IR properties of a 
sample of Seyfert galaxies as obtained with high spatial 
resolution images and found that the nuclear properties 
of the 12/im/[OIII] ratio showed no distinction between 
sources with or without a HBLR. This seems to suggest 
that despite the presence or absence of the HBLR, hot 
dust is still present in the nuclear region of both types of 
sources. 

Besides, Tran (2003) finds that although the IRAS 
25^m/60/^m color and the luminosity of the AGN are 
well correlated with the presence of a HBLR, the level of 
extinction towards the nuclei, is similar in both types of 
sources. 

We find some indications that sources with and with- 
out a HBLR might have systematic differences in their 
infrared emission. Bottom panels in Fig. |6] suggest that 
for sources with a HBLR the torus might be less extended 
than for sources with an undetected HBLR. However, as 
already discussed in § 5.2, the parameter Y is not well 
constrained by the fitting process of our data. 

Sources with a HBLR also might have systematically 
smaller line-of-sight extinction values than sources with 
an undetected HBLR. Eq. 2 states that the number of 
clouds along the line-of-sigh depends on the number of 
equatorial clouds (Ao), the thickness of the torus (er) 
and the inclination angle (Z.i). Fig. |6] shows that a is 
very similar for sources with and without a HBLR, while 
ZA does not show a statistically significant difference. 
The number of clouds, on the other hand, presents a 
more clear differences between sources with and without 
a HBLR. This can also be seen in Fig. [8] where nuclei 
with a HBLR cover a wide range of number of clouds 
at the torus equator, while nuclei without a HBLR clus- 
ter at the higher end of the distribution. However, Fig. [8] 
also shows that these differences are much less clear when 
looking at the number of clouds along the line-of-sight, 
A/"(/3), a parameter much closely related to the line-of- 
sight Ay (Eq. 3). Kolmogorov-Smirnov (KS) tests in 
fact do not confirm that the perceived differences are 
statistically significant. 

If further evidence that the different distributions of 
the line-of-sight Ay for sources with and without a 
HBLR are different is found, we need to explore some 
possible explanations. At face value this contradicts a 



strict unification scheme, where the only difference be- 
tween Type I and Type II sources is the view-point of 
the observer. However, is becoming clear that a strict 
unification scheme is not plausible (Elitzur 2012). 

We can postulate that a larger number of clouds ob- 
scures the polarized emission from the BLR. This would 
require a rather compact scattering region, with a compa- 
rable scale height as that of the torus itself. This result is 
in agreement with the analysis presented in Lumsden et 
al. (2004) where it is claimed that the fraction of Seyfert 
galaxies with a HBLR is larger when looking at only 
Compton-thin nuclei (as determined by X-ray observa- 
tions), implying a sample of sources with less obscured 
central regions. 

5.3. Torus Sizes and Masses 

The CLUMPY modeling does not provide an absolute 
torus size, but instead, the parameter Y = R ut/Rin- 
While most of the fits favor Y < 50, some sources can 
have very extended torus with Y < 100 (see Fig. [6|). 

Clumpy torus models are characterized by clouds show- 
ing a range of temperatures for a given distance from the 
central source (Nenkova et al., 2008b; Schartmann et al., 
2005). This is in sharp contrast with the predictions from 
a continuous dust distributions, where a unique temper- 
ature is found as a function of radius. Hence, the SED 
shapes for CLUMPY models are not very sensitive to 
the Y parameter or the size of the torus. Nenkova et 
al. (2008b) show that observations below 5 /im will not 
be able to distinguish between Y values, irrespective of 
the cloud distribution, which is determined by the q pa- 
rameter. At longer wavelengths some differences can be 
appreciated for a flat cloud distribution (q ~ 1), but 
these are only significant for wavelengths above 15/im 
and therefore are not well probed by our observational 
SEDs (see also §5.5). 

Following the results from Suganuma et al. (2006) we 
can assume that i?„ is indeed set by the dust sublimation 
radius as R ln = 0.4v/(X bol /10 45 ) pc, where L ho1 is in 
units of ergs/s. As already explained, we have used the 
[OIII] luminosities listed in Table |] as a prior for the 
intrinsic nuclear luminosities in our sources and found 
that the torus inner radii vary between 0.05 and 1 pc. 

Torus outer sizes (R ou t) are presented in Fig. [9] As it 
can be seen there is a wide distribution of R ou t , but most 
of them are below 5 pc in extent, with some torus being 
as small as 0.1 pc. 

The total mass of the torus can be estimated as 
(jNenkova et al.ll2008a[ ): 

M torus = Airm H sm(o-)N%R 2 in YI q (Y) (4) 

where is the equatorial column density of the 

torus, Ri n is calculated as the dust sublimation radius, a 
and Y are parameters of the model, and I q is 1 if q = 2 
or 3, Y/2lnY if q=l, or Y/3 if q = 0. We find some torus 
masses, up to 10 7 M Q , as shown in Fig. |9] However, in 
most cases these are driven by large Y values, a parame- 
ter not well constrained by the best fitting results. Still, 
Siebenmorgen et al., (2005) reported the dust masses im- 
plied by a simple model of the dust emission in 2 quasars, 
being of the order of 10 6 ~ 7 M Q . Fritz et al., (2006) re- 
ported on the masses implied by the smooth modeling 
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The median distribution is shown in solid color while the mode 
distribution is shown with a dashed line. 

of the emission of Type I and Type II Seyferts, ranging 
from ~ 70 — 1O 7 M . Our results show a strong peak, 
however, with most masses in the > 10 4 Mq range, as 
also seen by AH11. 

5.4. Comparison with Interferometric Observations. 

Infrared interferometric observations are providing the 
ultimate way to observe the putative dusty torus in AGN. 
Unfortunately, this technique is limited to the brightest 
sources and the observations required are very expensive 
and difficult to obtain. 

Until now long-based interferometry has been carried 
out for a couple of dozen Type I and II AGN (Wittkowski 
et al, 1998; Weinberger et al, 1999; Swain et al., 2003; 
Weigelt et al., 2004; Jaffe et al. 2004; Wittkowski et al., 
2004; Meisenheimer et al., 2007; Tristram et al, 2007; 
Beckert et al., 2008; Kishimoto et al., 2009ab; Raban et 
al., 2009), with the largest sample found in Tristram et 
al. (2009). Four objects in this last work are in common 
with our sample: NGC1068, F05189-2524, NGC5506 and 
NGC7582, however the observations for NGC5506 and 
NGC7582 did not provide useful data, while for F05189- 
2524 only a very faint fringe detection was possible. Also, 
it has been realized that interferometric studies yield 
more unambiguous results in Type I sources (see discus- 
sion in Kishimoto et al. 2011). Therefore, a one to one 
target comparison is not possible. 
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Fig. 10. — Radio- loudness, as defined by the parameter Rl ver- 
sus bolometric luminosity. The presence of a Hidden Broad Line 
Region is also indicated. 

However, for most of those sources where interfero- 
metric observations have provided restricting results, the 
sizes of the resolved structures observed in the near and 
mid-IR arc of the order of a few parsecs (see references 
above). These are already in good agreement with the 
results from the dust reverberation determined by Sug- 
anuma et al. (2006), where the inner face of the torus is 
found at the dust sublimation radius which directly de- 
pends on the luminosity of the central source, and with 
the results derived in this work. 

5.5. Radio-loudness 

Ho (1999) was the first to notice that radio-loudness 
is a function of the AGN bolometric luminosity, with 
low-luminosity objects (those below ~ 10 43 ergs/s) be- 
ing more likely to have a radio-loud central source. We 
can see this trend in Fig. [TUJ which includes all sources 
with L bo1 and Rl measurements. It can be seen that the 
probability of being radio- loud increases with L bo1 . 

It is interesting to see that for those radio-loud sources 
with acceptable fits (F 04385-0828 and NGC 7496) a very 
small inclination angle to the line of sight is derived from 
the SED results. This is in line with an orientation effect 
to be responsible for the boosting of the radio emission, 
although theoretical predictions state that the scaling be- 
tween the radio and the optical output sh ould depend 
weak ly on the relativistic Doppler factor (|Falcke et al.l 
2004). Further data would be needed to confirm this 
finding. 

In Fig. [10] we also include information about objects 
with HBLR detections. There seems to be a trend for 
more luminous nuclei to show the presence of polarized 
BLRs as already noticed by other works and in line with 
the hypothesis that weaker nuclear sources are out-shined 
by the stellar components (see references in §5.2). 

5.6. Correlation with Starformation 

One important subject to discuss is whether the pres- 
ence of nuclear obscuration is related to the level of star 
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Adequate fits were reported for 27 sources. 
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Fig. 11. — Luminosity of the 6.2/mi PAH feature, as a tracer of 
star forming activity, versus total obscuration at t he equator of the 
modeled torus. Symbols are the same as Fig. 1101 



formation in the nucleus vicinity. Taken at face value 
this scenario opposes the Unified Model, since in prin- 
ciple the only difference between a Type I and Type II 
source is the viewing angle. However, this is true only if 
all tori are exactly the same. 

As already discussed, Elitzur (2012) has recently pro- 
posed that Type II objects are examples of particularly 
puffcd-up tori, with larger A/o and a values. In turn this 
could result in a larger mass of cold dust, because of 
larger shadowing from the central source which allows 
for dust to cool more efficiently at large radii. 

In Paper I we have used diagnostic diagrams to deter- 
mine the presence of star-formation in the nuclear region. 
Essentially, there is no correlation between the level of 
starformation and the parameters that control the thick- 
ness of the torus, A/o and a. 

We can also look at the more extended star-formation 
using the observed luminosity of the 6.2/im PAH fea- 
ture. No evidence for such correlation is found between 
the strength of the PAH feature and the amount of ex- 
tinction determined from the SED fitting, as can be seen 
in Fig. HU 

6. SUMMARY AND CONCLUSIONS 

We have performed the fitting of the near and mid- 
IR SEDs of a sample of 36 Seyfert II galaxies using 
CLUMPY and 2pC models developed by Nenkova et 
al. (2002, 2008ab, 2010) and Stalevski et al. (2012). 



Our conclusions are the following: 

• Observations in the mid-IR, and in particular N- 
band spectroscopy of the sources, are crucial to 
perform an adequate fitting to the SEDs. 

• The use of the bolometric luminosity of the source 
as a prior during the fitting is also an important 
tool to constrain the best-fit results. 

• We find that the best-fit parameters for the 
CLUMPY models correspond to A/o > 5, a > 40, 
r < 25, Zi > 40, Y < 50 and A l ° B ~ 100 - 300. 
These values translate into typical torus sizes and 
masses of 0.1 - 5.0 pc and 10 4 ~ 6 M Q . 

• About half of the objects require an additional 
black-body component in the near-IR range to pro- 
vide an adequate fit to the SEDs. Most of best fit- 
ted temperatures are very high (T ~ 1700 — 2500 
K) and cannot correspond to the emission of hot 
dust. 

• 2pC models can sometimes provide a better fit to 
sources that require additional near-IR flux. 

• We find weak evidence that nuclei with HBLRs 
present lower levels of extinction than those with- 
out a HBLR. 

• Tentatively, we find that radio-loud objects are 
those with a very small inclination angle with re- 
spect to the line of sight. 

• We find no correlation between the torus properties 
and the presence of star-formation. 
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